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ie 6 SENTRODUCTION 


Many important structural problems in aerospace engi- 
neering involve components that are in contact with a fluid. 
One example is the fuel tank located in a wing or fuselage 
of an aircraft. Another example is the fuel tank in a 
liquid fueled rocket. An example in the field of civil 
engineering is the dam and reservoir systen. 

An ongoing research program at the Naval Postgraduate 
School is concerned with aircraft vulnerability in the fuel 
tank area. This program is currently devoted to a study of 
the Hydraulic Ram Problem. "Hydraulic Ram" is the dynamic 
loading of fuel tanks when impacted by bullets or warhead 
fragments. This loading gives the appearance of a simple 
pressure pulse with a resulting deformation or rupture of 
the tank. It has been found, however, that the loading and 
deformation are a combination of a number of events called 
"hydraulic ram components." [{Ref. 1] 

As part of the theoretical analysis of the fuel tank 
response to a penetrating projectile two computer codes were 
written to solve for the response of the fluid and the walls 
of a idealized two-dimensional fuel tank. The codes will 
compute the pressure throughout the fluid and will calculate 
the wall response to this pressure. The uSer can prescribe 


Eiewrontem tO be analyzed, including, but not 1AM eds 10: 





Fluid surface pressure disturbances, initial wall displace- 
nemo we Noid ewe vetOCt ples, pressure distributions through 
Tie mielital ae tex 

This thesis presents the method selected for the solution 
of problems involving two-dimensional fluid-structure inter- 
action. The types of problems which can be solved using the 
codes are described. Section II contains the theoretical 
description of the problem, and the following two sections 
outline two finite difference formulations that were used 
to obtain solutions. Section V contains a description of 
each code and describes the modifications necessary to 
change the boundary and initial conditions. Section VI 
presents the solutions obtained with the codes, including 
the check cases and several example programs, and discusses 
the comparison of results obtained from the two code formula- 
prons, Users Instructions and program listings for beth 


formulations are included in Appendices A and B. 





II. DERIVATION OF GOVERNING EQUATIONS 


A. THE WAVE EQUATION FOR AN INVISCID FLUID 


1. Equation of Motion 
Euler's equation of motion can be given in the form 


(i) 


+ grad p = O 


D 
oe 


where P is the fluid density, p is the fluid hydrostatic 


pressure (positive in compression), q is the velocity vector 
po is given by 


and t is time. The total derivative DE 


(2) 


Ug oe oO ms 
ue aes cme ay + w ae 


Dt ot 


where u, v, w are the components of the velocity in the 
Xx, yY, z directions respectively. When the fluid velocity is 


sufficiently small, 


(3) 


pm. 


=e 


She 
ct 


Superscript dots will be used to denote partial derivatives 


Wim respect tO time. 


Cons rucuseve Bquation 


(Ze 
The constitutive equation of the inviscid fluid is 


given by 
a= ke iy acl (4) 


where K is the bulk modulus. 





3. Equation of Continuity 
i Dp me = 3 
— ee 0 (5) 


Applying Eq. (3) we obtain 


ont 
He 
ae 


de GO (6) 
Eliminating div q between Eqs. (4) and (6) gives 


i (7) 


lick 
i 
DIP 
1S. 


4H. Potential Equation 


Define a potential 
q = grad 6 (8) 


Thus Eq. (1) becomes 


op grad 6 + grad p = 0 (9) 


when Eq. (3) is applied. Assuming p to be essentially 
constant throughout the fluid, Eq. (9) becomes 


grad (p¢ +p) = 0 (10) 


Ors 


od Fp = constant Gae 


0 





The constant term is set to zero since it.implies a base 
pressure which is constant throughout the domain. 


Differentiating Eq. (1) with respect to time gives 
pq + gradp = 0 (12) 


Applying Eqs. (4) and (8) 
5 grad $ + grad (-K div grad 6) = 0 @a) 


Since p is assumed to be essentially constant, Eq. (13) 


becomes 

erad(p¢ -~ K div grad ¢) = 0 (ie) 
or 

erad(po ~ KV~6) = 0 (15) 
Then 

a - KV" 6 = constant (16) 
or 

so > 
pd = KV oO (17) 


Recalling that the speed of sound in a fluid can be deter- 


mined from the relationship 


2_K g 
C ; Cie) 


albell 





_— 


we obtain the wave equation 


eek at = $ (19) 


B. THE GOVERNING EQUATION FOR THE STRUCTURE 
1. Two-Dimensional IJdealization 

Consider a long open box composed of five uniform 
plates. It can be idealized in two dimensions as a structure 
composed of three unit width beams attached at right angles 
as shown in Fig. 1. In the analysis each beam is considered 
separately, with its own set of independent boundary conditions. 
In all of the examples considered, each beam was assumed to 
be clamped at both ends.°™ 

The equation for the lateral displacement of 2 uniferm 


beam under dynamic load is [Ref. 2] 


w(x,t) + aie, =-) C20) 


tt] 
bay 
aq }& 
| 


where E is Young's modulus, I is the section moment of 
inertia, m is the mass per unit length, w(x,t) is the lateral 
displacement of the beam, and p is the pressure applied to 
the beam. A positive pressure is in the opposite direction 


as w(x,t) as shown in Fig. 1 for the tank. 


Ce) QNTERPACE CONDITIONS 
From Eq. (11) the pressure on the beam is the fluid 


pressure given by 


i 2 


HW a 
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pes = pd (238) 


where again a positive pressure is in the opposite direction 
as w(x,t). 

Another interface condition is the requirement that the 
normal velocity of the fluid potential at the interface is 


equal to the velocity of the wall, i.e. 


w(x,t) = $¢ (22) 
lal 


oP 


where a indicates a partial derivative normal to the wall. 
on 


Where the fluid has a free surface and the pressure is 


zero, the boundary condition is given by 
6 = 0 \23)) 


When an excitation pressure is applied to the free surface 


es p 
= fre Se 24 
$ ; (24) 
D. SUMMARY 
In summary, the equations used in the analysis are 
o“V"o = (19) 


14 
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III. NUMERICAL METHOD OF ANALYSIS — ACCELERATION FORMULATION 


A. FINITE DIFFERENCING SCHEME 
i. The Fluid 
The wave equation is hyperbolic in nature and for 
this reason an explicit formulation for the timewise response 
iemeciaescen. Ihe conventional central finite difference 
equations are used in both Space and time, and when applied 
to Eq. (19) give [Ref. 3] 


2 
Cc 


(Ax)* {645-1 7 os gt 7 464 3 y 54,3 i = 


ie 
At? {64-3 wie © bear} a (25) 


where Ax = Ay in a square mesh, i and j are mesh point 


reference indices, and k is the time step. is the 


os 54kt1 
only unknown. All other terms are known or can be evaluated 
prumber= te the caliculation of os 5 ,kt1 
2. The Beam Equation 
Conventional central finite dilference schemes are 


also used for the walls in both space and time. The resulting 


equation is given by [Ref. 3] 


16 





=. {w4_2 - Uw, , + 6, - 4wy,, + "s+2} ; 


= = + w = ~ 


where 
uae 
Paik 7 7 DAE {4.03 i" bear} | Sai), 


and j = 1lorje j 


A similar pair of equations can be written for the wall 


along the boundary where i= 1. In Eq. (26) w is the 


alee at 
only unknown in the explicit scheme. Note that the pressure 
Payee by FQ. (2/7) is a function of the fluid potential at 
Enemwall.. 
3. The Calculation Scheme 
When the two systems are combined, the other linking 
mechanism is the interface condition on the fluid and wall 


velocities. Using central finite differences, Eq. (22) is 


eivem by [Ref. 3] 


1 1 
Fare — = —— om + 
2At Wiyker ” "5 cer} 2hx oS cail gl ‘_e ai 


Ly 





In the calculation scheme the wall deflection at each new 
step ktl is computed first, and then the velocity potential 
MmomCOMmpuLea tnrouscnouct the fluid. In order to:compute the 
wall displacement at k+l the pressure must be known at the 
preceding time step k, which requires knowledge of the 
velocity potential at the wall at the new time step, ktl, 
according to Eq. (27). The solution to this problem is to 
solve for the wall displacement by implicitly solving for 
the new velocity potential at the wall. Thus, using 


Eq. (25) to eliminate 9. from Eq. (27), and using 


Ls deaketel 
Eq. (28) to eliminate d. cot ie (when j=1) from the result 
2 > 


gives 





+ _ . + ae 
Sean see 4k * PL j+1sk 


Ax 
* Or jti,k * At (Ws J} 


A similar equation can be obtained for the wall where 
j= j MES TOS tTencine pesiven by EG. (29) into Eq. (26), 


combining terms, and solving for w. gives 
eet 4: 
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- Pp er 
Wa k+l bm {264 5 eos jyk-1 7 


cAT\,2 
caw ae = 
ie Mea ttt 7 Papigie t. 2°4 pox 


Ax 
re a); 





Week ~ 4 K-21 
a At ie W - wy + 6w W 
mM *Cax)e TD te ibe i,k it1,k 
+ w e+ i até a) 
i4+2,k m Ax (30) 
Equation (30) is written for the wall where j = 1, 0 < x < aH. 


For x > aH, p = 0. Once the wall deflection is known at ktl, 


the interface conditions on the fluid are known. The calcula- 


tion of the fluid velocity potential is now quite straight- 
forward using Eq. (25) and Eq. (28) at an interface. 

On the free surface of the fluid the velocity 
potential may be excited by an applied pressure. In this 


case, applying central finite differences to Eq. (24) gives 


ies At 


where (Py), is the pressure exerted on the boundary at time 


step k. 


Ie 








B. INITIAL CONDITIONS 
1. The Fluid 
Since the governing partial differential equation 
(Eq. (19)) is second order with respect to time, two initial 
conditions are required; 1) the velocity potential at time 
zero, and 2) the first time derivative of the velocity 
potential, which is Sropenttona) to any initial pressure in 
179 0X a bo 
2. The Wall 
The beam equation (Eq. (20)) is also second order 
with respect to time, and two initial conditions are required; 
inmemc Initial Gispllacement of the wall, and 2) the initvwal 


velocity of the wall. 


eee the ive ore POR NUMERICAL STABLLITY 
1. The siuid 
For hyperbolic systems solved using Eq. (25), the 


non-dimensional time step ratio 


aged 
Eee (32) 
where primes denote non-dimensional quantities, must be less 
than 1 for numerical stability for a one-dimensional system 
[Ref. 3] and must be less than {2/2 for numerical 


Stability for a two=diménsional system [Ref. 4]. 
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Non-dimensional time t' is given by 
re LS 
t 7 (33) 


where L is a reference length. The non-dimensional length 


x' is given by 
xt = & CER), 


Substituting Eqs. (33) and (34) into Eq. (32) gives 


= Ax 


OX (35) 
oun oe 


for numerical stability of the two-dimensional system. 
2. The Wall 
Non-dimensionalizing the time and length coordinates 


in the beam equation gives [Ref. 3] 


ea: Pee 
t' = = a (36) 
x' = : (37) 


The beam equation is parabolic in nature [Ref. 3] and 


requires a non-dimensional time step ratio 


__ Att ol 
* Caxty? = 2 oy 


eu: 





for numerical stability. Substituting Eqs. (36) and (37) 


into Eq. (38) gives 


ue Dac bi 
So ae) es (39) 


Therefore 





- 
ar < Me (40) 


3. The Combined System 
The assumption has been made that the fluid will 
require the smallest time step for stability. Thus, there 
is a restriction on the value of the parameters used in the 
combined analysis, in particular the physical parameters of 


the beam. This restriction is given by 


a Ax ee _ 
es mle Chale) 
fo —- 2 Ba 


according to Eqs. (35) and (HO). 
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IV. NUMERICAL METHOD OF ANALYSIS — VELOCITY FORMULATION 


A. VELOCITY FORMULATION APPROACH 

Kreig and Monteith [Ref. 5] have suggested a modification 
to the timewise finite differencing Scheme outlined in 
Section ©l we elhis Modimicarton 1s Cescrivedecoucus Ve loci ia 
POremulct ON. 

Consider their example of the one-dimensional spring mass 


System. 
Mw + Kw = 0 (12) 


Using central finite differences, Eq. (42) becomes 





M { } 5 
W -~ Pw. + Ww Fea = 0 (43) 
(At)= k-1 Fe k+1 k 
or 
B: K 2 
Now consider the velocity formulation 
My + Kw = 0 
w oF Vv (45) 
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Defining the velocity at half time steps and the displacement 


at whole time steps leads to 


= {- ver + Victig} ate KW), = 0 (46) 
or 

“KHs VU k-3s 7 Ot my, (47) 
and 

Seg, Me eNO Who (48) 


With either formulation the time step is the same, and in 
the real number system without round-off error the results 
would be identical. However, computers do not compute in 
the real number system. Round-off errors are always present. 


Consider the term in brackets in Eq. (44) 
K o 
[2 - M Gite) = || 


Depending on the number of digits retained, the significance 


of = (At) could be completely lost in round off. 


au 





beet INITE DIFFERENCING SCHEME 
1. The Fluid 
The fluid equation derived in Section III can be 


written as 
c“V%o = (49) 
$= (50) 


A central finite difference scheme is proposed centering 
the "velocity" at whole time steps and the "displacement" 
at half time steps. Note, the pressure is directly computed 


at each whole time step since 
Py; a eas ©, Vy, (ay) 


The finite difference equation for the fluid is given by 


i {Veo : vic} _ ae 


where ys ; oy is the only unknown. The potential is obtained 
= 3 2 


HeOMmuner minite difference expression for hq. (50) 


+ Atw, (53) 


Os j,kts ~ 94,5, ks j5k 
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2. The Wall 


Similarly, the beam equation can be written 


y 
lal an w(x,t) + mv = -p (54) 
ox 
and 
w(x,t) =v (55) 


For these equations a central finite difference scheme is 
proposed centering the velocity v at nalf time steps and 
the displacement w at whole time steps. Note that the 
velocity needed for the interface condition is directly 


computed, 


= 2s — —“s 
“i,k-s Ox ae assaf}, ; (56) 
aed | 


The finite difference equation for the beam is 


m = 
At a : Vi tl P Ws 54k ad 
The unknown is Va ths From Eq. (55) 
Ms yet = “ak + At Vath (58) 
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peaeene Combined System 
With the velocity function there is no need to solve 
for the beam displacement in terms of an implicit solution 
of the potential at the wall. In the velocity formulation 
the pressure is computed directly; no differencing is 
required. Figure 2 graphically illustrates the time step 
centering of both the acceleration formulation and the 
velocity formulation. 
The solution scheme for the velocity formulation is: 
1. Solve for the wall velocity, Eq. (57), based on the 
known pressure at the wall and the previous displace- 
ment and velocity. 
2. Calculate the new wall displacement, Eq. (58), based 
on the new wall velocity and the previous displacement. 
3. Solve for w, Eq. (52), based on the previous wy, the 
new wall velocity (for the interface conditions) and 
the previous 6. 
4, Calculate the new $, Eq. (53), based on the previous 
@ and the newly calculated vy. 
See Return to ae 1 and repeat. 
As in the example given by Kreig and Monteith this formula- 
tion is exactly equal to the acceleration formulation if the 
beam and fluid equations are considered separately. However, 
in the combined system the fluid equation is displaced 
exactly one-half time step from the fluid equation used in 


The acceleration formulation. In the acceleration formulation 


aay 





ACCELERATION FORMULATION 


QE 


Oo = 


TIME 


a) 


+2 


+1 


VELOCITY FORMULATION 


fx] 
S 
i 
FA 
———_, 
pss 
— 
| a 
11 = 
| ot ee 
|) > | 
pes | | 
Nay ] 
t+ ce ] 
Pd 
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ase > | 
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ee el ee 
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eal 


de GY 
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mae pressure depends on Ort] and dy-1 Whereas iimeume velocity 
formulation, it depends on ee and pee ». Thus, even with 
no round-off error there will be a difference between the two 
solutions. Refer to Part D for a discussion of the signif- 


ieanece of this difference. 


C. SNITIAL CONDITIONS 

As with the acceleration formulation a total of four 
initial conditions are required. deere recall that two 
of the four variables are defined one-half time step on 
either side of t = 0. They are the wall velocity and the 
Woloec uy pevencial Of @ In uhe fiuad. 


The actual initial conditions can be denoted as Wo and 


dé, : The computer code assumes that 
Vv = 5 
wy, ee (59) 
and 
= 60 


Bs JHE TIME STEP FOR NUMERICAL STABILITY 

Theoretically, since the same parameters are used in 
both formulations, the time step for numerical stability 
Should be exactly the same for both formulations. In fact, 
this was found to be true for the beam equation and for the 
fluid equation when each was tested Separately. However, 


Maen the Systems were coupled, and.the maximum theoretical 


ao 





time step was used, substantial instability exhibited 
itself immediately in the computer results. The instability 
occurred at the interface and it occurred within two time 
steps of a pressure pulse striking awall. A subsequent 
velocity was imparted to the wall which was always large 
enough to cause a reversal in the fluid pressure at the wall 
at the next time step. 

A primitive but effective solution to the instability 
problem was obtained by reducing the time step arbitrarily 


Ee 
Ne @ es (61) 


0 13 
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V. DESCRIPTION OF THE COMPUTER PROGRAMS 


Two Fortran IV programs were written and tested. The 
theory behind each is outlined in Sections III and IV above. 
The two ante ae are quite similar in construction and 
require similar amounts of core storage. C.P.U. time 
required is also similar. Only one program will be described 
in detail, but differences between the two will be pointed 


Out where they occur. 


A. GENERAL FEATURES 

The programs allow for complete flexibility in any of 
the possible initial conditions, both on the walls and on 
the fluid. At present the program ailows for three walls 
cea Varaaple fluid Gepin, @ihe ends Of each wall. are 
considered to be clamped. The wall boundary conditions can 
be changed by any future user by changing six equations in 
the velocity formulation and twelve in the acceleration 
formulation. The location of these equations will be pointed 
out below. The program provides for a non-dimensional thirty 
by thirty mesh point domain, with the reference length being 
the length of the horizontal wall (See Fig. 1). All variables, 
vectors and arrays are placed in blank Common, to allow for 
complete flexibility in communication between subroutines. 
Purcnnew all vectors, and arrays are ‘hagueaansuped vorZerO Value 


before analysis begins. 
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Subroutines are provided to: 

1. Compute constants derived from input parameters, 
normalize the dimensions of the problem, and compute 
or sev initial @eondivr1 ons. 

2. Output data obtained for each mesh point in the domain. 
The format is tabular. 

3. Compute the beam equations. 

1. “Compute thestluidueguetions. 


5. Provide for a special graphical output. 


Each subroutine will be described in detail. 


B. MAIN PROGRAM 


imeevorlable Desemiptiens sv clocsry=e! ermular1 On 
WA(T) - left vertical wall displacement at k 
WB(T) - horizontal wall displacement at k 
WC(T) - right vertical wall displacement at k 
VA(T) - left vertical wall displacement at ktl 
VB(T) - horizontal wall displacement at k+l 
iC Cis) - right vertical wall displacement at k+l 
UA(T) - left vertical wall velocity at k-'s 
UB(T) - horizontal wall velocity at k- 
Wiens) - right vertical wall velocity at k- 
XACL) —- left vertical wall velocity at kt 
pds (OL) - horizontal wall velocity at kt 
xCGE) - right vertical wall velocity at kt 
WADOT (I) - initial left vertical wall velocity 
WBDOT (T) - initial horizontal wall velocity 





WCDOT(TI) 
PRA(T) 
PRB(TL) 
mic L) 
PK(T) 
IWA (I) 


IWB(T) 


IWC (TL) 


P(E sd) 
PN(I,J) 
Peeler) 
PRESS(I,J) 
PHIDOT (I,J) 


IPRESS(I,J) 


ier right vertical wall velocity 


initial left vertical wall pressure 


initial horizontal wall pressure 


initial right vertical wall pressure 


fluid free surface applied pressure 


integer value of left vertical wall 
displacement 


integer value 
displacement 


integer value 
displacement 


value of ¢ in 


value of ¢@ in 


value of win 


of herigontel wall 


of right vertical wall 


the -tluld at k= 
the fluid at kts 


the fluid at k-l 


value of pressure in fluid at k 


initial value 


integer value 


In the displacement formulation 


VA,VB,VC 


WA,WB,WC 
dee 


UA,UB,UC 


P(I,J) 
PN(I,J) 


EG) 


are not used 


au G 


ace 


are 


aie 


are 


aice 


values 
values 


values 


values 
values 


values 


Os 
of 


of 


of 
Os 


One 
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CU in the fiuid 
of pressure in the fluid. 


the following variables differ. 


wall displacement at k 
wall displacement at ktl 


wall displacement at k-l 


6 at k 
dat ktl 


@ at k-l 





2. Main Program Functions 
In the velocity formulation the main program 
accomplishes the following functions: 
Mm Conurols Calls To Ouhner Subroutines. 
ee Computes first value of PN(IlJ) from initial conditions. 
3. Saves the value of the displacement at a selected 
wall mesh point at selected time steps for later 
plotting as a time vs. displacement raphe. 
4. Turns off a free surface pressure disturbance at any 
preselected time step. 
In addition the "acceleration" formulation program 
5. Computes the k = 1 value of wall displacement from 


iMieaial COndLtions. 


C. SUBROUTINE INP(IER) 
this subroutine 
ieee Reads ail input data. 
EeeeecOmpuLes constants. 
Pee cote vatues TO initial conditions other than zero 
Wome: 


i. Variable Description 


RHOF - Density of the fluid 

E - Young's Modulus 

RL ~ Horizontal length 

RH sevemoreal height where RH < RL 

AL - Fractional depth of fluid 0 < AL < 1.0. 
C ~ Speed of sound in the fluid. 


34 





did at 
RHO 


NH1 


N1 


DELX 
DELT 


RIO 


NA 


e15C2,C3 


ee). Co,c;,Co,co,D 


D1,D2,D3 
NH 


2, Inpueecards 


Thickness of beam 


Densasty OClapeam material in ibs 
per in 


Number of vertical mesh points 


Number of horizontal mesh points 
Wieiecn Nits se Ni = 730 


Width of mesh 
Time step 


Computed value of I, moment of 
inertia 


Number of vertical mesh points in 
thet turd. 


ConsGants -cemmon to both programs 


Constants used only in acceleration 
ToOrmibar Lon 


Constants used in velocity formulation 
Equals NH1 - l 
Equals NH1 - 2 
Equals Nl - 1 


Equals NA - 1 


OniveLWO InDUGEcards are tised vO input Rie. RH, AL: 


C, TH, RHO, RHOF, E. 


Appendix A. 


The format will be described in 


3.  Non-Zero Initial Conditions 


Any non-zero initial conditions are inserted by the 


WSer "ac the end of Subroutine INP. 





D. SUBROUTINE OUTP(IER) 
Subroutine" OUTP™ prints out in tabular’ format the values 
of each wall mesh point displacement and the pressure at 


each mesh point throughout the domain. 


E. SUBROUTINE WALL( IER) 

Subroutine WALL computes the values of the wall velocity 
and displacement in the velocity formulation and the wall 
displacement in the acceleration formulation. If a future 
user desired to change the boundary conditions of the beam 
equation, then any equation with a computed answer stored 
in a variable indexed with a 2, N, or NH would have to be 
changed. In the velocity formulation only the equations 
for VA(2), VA(NH), VB(2), VB(N), VC(2), VC(NH) need to 
be changed. In the acceleration formulation these equations 
are for XA(2), XA(NH), XB(2), XB(N), XC(2), XC(NH) and recall 


that these equations must also be changed in the main program. 


Eo UCROULINE FLUIDC IER) 

Subroutine FLUID computes the value of YW and ¢ at k and 
k+¥5 respectively in the velocity formulation and computes 
Gmeevalues of 6 at k+l in the acceleration formulation. 


Eregsmre at each mesh point is also computed in FLUID. 


G. SUBROUTINE GRAPHO(IER) 
Subroutine GRAPHO is an attempt to make some Sense out 
SumunesmaoOssible 1900 numerical values calculated at each 


time step. The subroutine locates an absolute maximum value 
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for wall displacement and fluid pressure and from these 
values normalizes the entire domain to integer values from 
-99 to +100. This information is then presented in a format 
quite similar to Fig. (1). An integer value is Placed at 


each mesh point in the domain. 





VI. SOLUTIONS 


A. CHECK CASES 

In order to determine if the programs correctly modeled 
the system described above, several test cases were run on 
each formulation. 

1. .Lheshdaad 

By controlling the subroutines called by the main 
program, the user is able selectively to compute results 
for the fluid alone or the wall alone. One test case for 
the fluid alone was a step pressure disturbance applied 
uniformly over the entire upper surface for a finite time 
interval. The walls were treated as rigid. This test 
showed that the pulse propagation and reflection behaved 
as expected in both formulations. The uniform pulse traveled 
uniformly through the fluid. The value of pressure doubled 
as expected at the horizontal rigid wall and returned 
uniformly to the free surface. See Figs. 3 and 4. 

Another test case consisting of a point pressure 
Gisturbance at the middle of the free surface Showed that 
two-dimensional pulse propagation in the fluid was as 
expected. The amplitude of the pressure decreased radially 
from tle peint source. 

In both test cases the theoretical time step for 


Stability proved wo be valid in both formulations. 
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2. “ihe Wall 

A free vibration analysis was conducted using the 
beam equation to test the wall response alone. The wall 
was initially displaced to conform to the first mode of 
vibration. Rogers [Ref. 1] gives a value of 51.3 cycles per 
second for the value of the frequency of the first mode of 
vibration for the beam parameters used in the test. The 
results from both formulations returned a first mode of vibra- 
tion frequency of approximately 50 cycles per second. Again 


the theoretical time step for stability proved to be valid. 


B. EXAMPLE PROBLEMS 

The types of problems these two codes were developed to 
solve involve the combined system. The first example prob- 
lem consisted of a uniform pressure pulse at the free sur- 
face of a flexible tank for 2.26 milleseconds. When a solu- 
tion was first attempted using the velocity formulation, the 
previously mentioned instability was discovered when the 
theoretical time step for numerical stability of the fluid 
was used. Since no other solutions to this type of problem 
were available, first priority centered on determining the 
Cause and cure of the instability. A reduction in the time 
step as discussed in Section IV solved the instability prob- 
lem. The velocity formulation was run with a reduced time 
Step, and the acceleration formulation was run using the 


theoretical time step. The results for the displacement 
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of the center of the horizontal wall are shown in Fig. (5). 
The dominant frequency of vibration is approximately 12 
cycles per second for both formulations. There is a small 
amplitude discrepency of about 5% between the two formula- 
tions. The resulting vibration appeared to involve all 
meades in the beam. 

The acceleration formulation was also used to solve the 
same problem with a uniform pulse held for .226 milleseconds 


or 1/10°% 


the previous example. Figure (5) again denotes 
a frequency of approximately 12 cycles per second, and the 
amplitude is proportional to the previous case. Maximum 
amplitude is about one-tenth that of the previous case. 

It is interesting to note that if the fluid is considered 
as incompressible, and its mass is added to the mass cf the 
beam, the frequency of free vibration as given by Rogers 
[Ref. 2] is 11.1 cycles per second. 

Another example problem was run with a pressure exerted 
at the two center mesh points of the free surface for approx- 
imately 450 microseconds. In this case, the first mode of 
vibration predominated in the horizontal wall at a frequency 
of 13 cycles per second. Again only minimal differences 


in amplitude were observed between the two formulations. 
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VII. CONCLUSIONS 


The results obtained to date appear to support the 
following views: 

1. The amplitude of vibration in the walls are depen- 
dent upon the length of time the free surface excitation 
pressure is applied. This would be expected purely from 
energy concepts. 

2. Choice of formulation is up to the user. Each has 
its own virtues and its vices. The velocity formulation 
requires less core storage and less CPU time for the same 
number of time steps. (The acceleration formulation requires 
15% more CPU time than the velocity formulation; however, 
the time step required is 19 percent smaller in the velocity 
Pormuacvion.) The velocity formulation would be easier to 
moar y. 

3. No claims are made as to the efficiency of either 
pregram as presently written. In fact it is estimated that 
the program size could be reduced considerably. 

4, More work needs to be done in checking out other 
simple cases to prove the validity of the codes, for exam- 
ple, pressure pulses passing through the fluid at angles 
Peonmune Horizontal. . 

Se otuener program will @ilve-a fusure uSer valuable 


insight about structural behavior in contact with a fluid. 
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APPENDIX A 


USER'S INSTRUCTIONS 


Input of initial conditions has already been discussed 
in Section V above. Actual input data is contained in two 


cards with the following format: 


Card 1 
Column Variable Definition Format 
Sara il: 
1-10 RL Reference length F10.4 
11-20 RH Vertical height Ome 
21-30 AL Fractional depth F10.4 
Olt toad 

31-40 C Vemloc. ty 10f ssound F10.4 
Ln cient UG 

41-50 | Thickness of wall F10.4 

51-60 RHO Density of wall 3 F10.4 
material in #/in 

Cara 2 

1-12 RHOF Density of fiuid liz 
in #-sec e/in 

13-24 E Young's Modulus of mice 


wall material 
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